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ABSTRACT 



In late 2006, ground-based photometry of 77 Car plus the Homunculus showed 
an unexpected decrease in its integrated apparent brightness, an apparent rever- 
sal of its long-term brightening. Subsequent HST/WFPC2 photometry of the 
central star in the near-UV showed that this was not a simple reversal. This 
multi-wavelength photometry did not support increased extinction by dust as 
the explanation for the decrease in brightness. A spectrum obtained with GMOS 
on the Gemini-South telescope, revealed subtle changes mid-way in rj Car's 5.5 yr 
spectroscopic cycle when compared with HST/STIS spectra at the same phase 
in the cycle. At mid-cycle the secondary star is 20-30 AU from the primary. 
We suggest that the spectroscopic changes are consistent with fluctuations in the 
density and velocity of the primary star's wind, unrelated to the 5.5 yr cycle 
but possibly related to its latitude-dependent morphology. We also discuss sub- 
tle effects that must be taken into account when comparing ground-based and 
HST/STIS spectra. 

Subject headings: General, Line: Profiles, Stars: Individual: Constellation Name: r] 
Carinae, Stars: Variables: Other, Stars: Winds, Outflows 
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Introduction 



Eta Carinae has a peculiar habit of settling into a photometric trend and spectroscopic 



state long enough for observers to grow accustomed to it, and t 
Dramatic transitions occurred around 1840, 1890, and 1945, see 



ren abruptly changing 



Humphreys et al 



teoosh 



and refs. therein; and there are hints that another may have begun in the late 1990's. For 



almost 50 years after the 1940-1950 c hange, the star p 



of 0.025 visual magnitude per year (lO'Connell 



1974 



Mattei & Foster 



1998 



Davidson et al. 



1956 



1999a 



us eject a 



Feinstein 



Martin 



jrightened at an average rat e 



967 



Feinstein fc Marracc 



20051 ). largely due to expansion 



of the dusty Homunculus nebula which had been ejected around 1840. Essentially a 
bipolar reflection or scattering nebula, the Homunculus dominated r\ Car's total brightness 
throughout the 20th century. After 1997, however, observations with the Hubble Space 



elescope (HST 



( Davidson et al. 



revea l ed that the central star was now brightening far more rapidly 



1999a ; 



Martin fc Koppelman 



2004 



Martin et al. 



2006b|). Despite brief 



declines during the 1998.0 and 2003.5 spectroscopic events, its average rate from 1998 
to 20 06 was about 0.15 mag/yr. This was a recent developm ent, since earlier speckle 



data ( jWeigelt fc Ebersberger 



1986 



Hofmann fc Weigelt 



19881 ) show that the star had not 



brightened that quickly in the fifteen years before 1997. Ground-based photometry of the 



Homunculus after 1997 showed an increase that was less dramatic, 



of the fluctuations seen between 1955 and 1995 (Sterken et al. 



Fernandez-Lajus et al 



2003 



Whitelock et al 



2004; 



j ut which exceeded any 



1999 



Davidson et al. 



Martin fc Koppelman 



1999a; 



2004T ) . Most 



likely the circumstellar extinction began to decrease rapidly in the mid- 1990 's; perhaps 
the rate of dust formation near the star had decreased, or dust was being de stroyed, or 



both. Either case requires some undiagnosed change in the dense stellar wind (IMartin et al. 
2006tj). 



Thus, suspecting that another "change of state" was underway, we were surprised 
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when ground-based photometry of rj Car plus the Homunculus showed a 0.3-magnitud e 



decrease in brightness 



Fernandez-Lajus et al 



2008 



between August 2006 and March 2007 ( DiScala fc Jones 
20091 ) J 1 ! The change occurred at all BVRI wavelengths and initially 



resembled the behavior seen during a spectroscopic event (Fig. [T]). This was the largest 
decline in r? Car's grqundb a sed light curve e vent in at least the past 50 years, see 



van Genderen et al. 



(1999) 



Davidson et al. 



(Il999af ). and refs. t herein. The integrated 



photo metry of the Homunculus and central star vary in parallel ( Ivan Genderen fc Sterken 



2004 ) . Therefore, if this change proved to be a reversal of the previous several years' 
behavior, then one might expect the central star to have faded by more than half a 
magnitude in 2006-2 007, largely invalidating our conjectures about the future trend 



( Martin et al. 



2006bf ) . Thus it was important to obtain fresh HST photometry of the central 
star. Fortunately this was already a goal of a small existing program outlined in Section 2 
below. 

As usual in rj Car research, the results differ from any scenario proposed beforehand. 
The photometric decline soon leveled off, but our HST photometry strongly suggests that 
it was not merely a temporary reversal of the previous brightening trend (§3). Suspecting 
that the photometric reversal might signal unforeseen spectroscopic changes in the stellar 
wind, we obtained data with the Gemini Multi-Object Spectrometer (GMOS) in 2007 which 
can be compared with earlier data from the HST's Space Telescope Imaging Spectrograph 
(STIS) in 1998-2004, and with VLT/UVES spectra from 2006 prior to the photometric 



See the CCD photometry by Georgio DiScala (Australia) and by Raymond W. Jones 



(South Africa) in the AAVSO International Database, http://www.aavso.org/ and the pho 



tometric record from La Plata Observatory (Argentina), http://etacar.fcaglp.unlp.edu.ar/ 



In our Figure [TJ, note that the decline may have begun in the season that is unfavorable for 
observing rj Car. It may have started as early as August or as late as December of 2006. 
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decline. The spectroscopy revealed only subtle changes halfway through rj Car's 5.5-year 
spectroscopic cycle, possibly related to the latitude-dependent morphology of the stellar 
wind. 

Any interpretation of i] Car's recent beha vior must take into ac c ount it s 5.54-year 



Humphreys fc Stanek! (120051 ) . 



Martin et al. 



spectro scopic cycle; for general information see 
( I2006a[ ). and many references cited therein. Most authors agree that the 1998.0, 2003.5, and 
2009.0 "spectroscopic events" probably occurred at periastron passages of a hypothetical 
companion star in a highly eccentric orbit. At times more than 7 months from such an 
event, the two stars are 15 to 30 AU apart - seemingly too distant for the secondary star 
to significantly influence the appearanc e of the primary wind in the observations discussed 



here. Although 



Kashi fc Soker 



( l2008af l have suggested that even at apastron the ionizing 
radiation from the secondary may not be negligible. At least in most scenarios proposed 
so far, one expects the 2006-2007 changes to signify structural effects in the primary wind 
more than wind-wind interactions. 

In the next section we describe the photometric and spectroscopic observations 
including a description of the differences between ground-based and HST/STIS spectra and 
the complications that must be accounted for when comparing them. In §3] we discuss the 
observed variations in n Car in 2006-2007, mid way in its 5.5 yr spectroscopic cycle, and 
compare them with spectra from the same phase in the previous cycle (2000-2002). The 
last section summarizes the observations with a brief discussion of their possible origin. 
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Observations 



2.1. HST Photometry 

Normal ground-based images of rj Car do not consistently separate the star from the 
inner parts of the Homunculus nebula, which has structure on all scales from 0.1" to 10". 
Until recently the Homunculus appeared brighter than the st ar, though the ejecta -to-star 



2006bh . Precise 



brightness ratio has evolved substantially in the past 20 years (IMartin et al.l 
ground-based photometry generally includes the entire Homunculus plus the star. HST 
images, however, have good enough spatial resolution to show the brightness of the star 
itself with only minor contamination by the ejecta located at r > 0.15". In recent years we 
have obtained a number of images of this type, to monitor the central star's brightness and 
detect changes in the surrounding ejecta. 

We obtained Hubble Space Telescope ACS/HRC images of r\ Car in the F250W 
and F330W filters on 2006 August 4 (MJD 53951.1), near peak brightness, and 2007 
January 20 (MJD 54120.4), during the unexpected decline. The central wavelengths of 
these broad-band filters are near 250 and 330 nm respectively. All of t he images were 



reduced, measured, and calibrated as described in 



Martin et al. 



(l2006bl ). For brightness 



measurements we used a weighted virtual aperture centered on the star, with diameter 0.3" 



(about 10.7 AC S/HRC pixe 







.15"S 



See also 



Smith et al. 



s ). Th e aperture weighting function was 1 — r 2 /R 2 with R 
(120041 ) for aperture size effects in ACS images of r\ Car. 



Shortly after our January 2007 observations the ACS electronics failed. We then used 
the older HST/WFPC2 camera on 2007 August 23 (MJD 54335.1) and 2008 February 14 
(MJD 54510.8), with filters F255W and F336W which are comparable to the ACS/HRC 



- 



Tf 



2R. 



In IMartin et al.l (I2006bf) we wrongly quoted R = 0.3", which was really the diameter 



Martin fc Koppelman 



(120041 ) stated this detail correctly. 
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filters. These images were reduced using the standard STScI data reduction pipeline. A 
correction factor of 0.9915 was app lied to the raw counts to compensate for geometric 



pixel distortion (IBiretta et al. 



19961 ). and the PHOTFLAM keyword induced a pipeline 



conversion from raw counts to flux values. 

We measured the brightness of the central star in the WFPC2 data, using the same 
0.3" weighted virtual aperture (~ 6.6 PC2 pixels) as for the ACS/HRC data. We estimated 
the aperture correction based on two 12th-magnitude stars in the images, well separated 
from the Homunculus, by comparing their fluxes in the virt ual aperture to those in a 



1995|). The 



standard 1" measurement aperture corrected to infinite size (IHoltzman et al. 
resulting geometrical transmission factors are listed in Table [TJ 

To assess the effect of a known "red leak" in the WFPC2 F336W filter transmission 
curve, we convolved an HST/STIS spectrum of rj Car with the response function for the 
Planetary Camera with this filter. The flux at A > 450 nm contributed only 0.94% of the 
counts, negligible for our purposes. 

The STMAG photometric system is cal ibrated for direct co mparison of fluxes measured 



with similar filters in different instruments ( Sirianni et al. 



20051 ). We tested the calibration 



by convolving the instrument and filter response functions with HST/STIS spectra of rj Car 
and comparing the results for the ACS/HRC F250W and F330W filters to those for the 
WFPC2 F255W and F336W filters. The computed values agreed to four significant digits, 
giving us confidence that comparisons between our ACS/HRC and WFPC2 results are 
valid. Our experience with the large existing HST data set for rj Car, 1991-2009, has never 
contradicted this assessment to a serious degree, despite the spectral forest of emission lines 
and strong UV absorption. 

The ACS/HRC and WFPC2 photome try are listed in Tab le |2l and plotted in Figured] 
together with our previous measurements (IMartin et al.ll2006bl ). 
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2.2. Spectroscopy 



We obtained ground-based slit spectroscopy of rj Car with the Gemini Multi-Object 
Spectrograph (GMOS) on the Gemini South telescope in late June and early July 2007, 
just after the decrease in brightness. The parameters of these observations and subsequent 
GMOS observations in February and July 2008 are listed in Table [3j The parameters of 
the spectrograph are given in Table HI We prepared 2-D spectrograms with the standard 
GMOS data reduction pipeline in the Gemini IRAF pa ckage and extracted 1-D spectra 



via a routine developed earlier for use with HST/STIS (jMartin et al 



2006af ) (Figure |2). 



At each wavelength our software integrates the counts along a line perpendicular to the 
dispersion, weighted by a mesa-shaped function centered on the local spectral trace. For 
these data we used a mesa function with base- width = 11 pixels and top- width = 7 pixels, 
about 0.8" and 0.5" respectively. The seeing was roughly 1.0", so each GMOS spectrum 
discussed in §3 below represents a region about 1" across. 



HST/STIS spectra of rj Car had been obtained in October 2001 and 



2002 , at approximately the same phase in its 5.54-year spectroscopic cycle (IDavidson 



anuary 



2004 ). We reduced t hem using a modified version of the Goddard CALSTIS reduction 



pipeline described in 



Martin et al. 



( 12006a] ) . The spectra are listed in Table [3] and the 



instrument parameters are given in Table HI The January 2002 data most closely match the 
spectroscopic phase of July 2007^ 

Of course the spatial resolution of Gemini, limited by atmospheric seeing, is greatly 
inferior to HST/STIS and the inner ejecta are unavoidably mixed with the spectrum of 
the star. In particular our GMOS spectra include the 'Weigelt knots" 0.3" northwest of 



3 STIS and GMOS spectra used in this paper can be downloaded at 
|http:/ /etacar.umn.edu/ 
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the star. Fortunately, by 2007 the star had be come substantially 



ejecta - unlike the case 10 or 20 years earlier (jMartin et al 



Drighter than the nearby 



2006ah . The 1-D STIS spectra 



used in this paper were prepared with an extraction width of 0.20", which excludes most of 
the narrow emission from the inner ejecta; see important comments in subsection 12.31 

Moreover, the slow-moving inner ejecta produce narrow emission lines which are 
distinguishable from the broad stellar wind lines; typical widths are of the order of 20 and 
400 km s~ x respectively. In the wavelength range of interest, the spectral resolution of 
GMOS is roughly 75 km s _1 while STIS provided about 40 km s -1 . The narrow lines are 
therefore more blurred in the Gemini data while the broad stellar wind features and their 
P-Cygni absorption components are well resolved by both instruments. We have smoothed 
the STIS data plotted in Section [3] to match the resolution of GMOS. 

Ha emission is so bright in rj Car that it saturates the detector pixels in even the 
shortest GMOS exposures centered on the star. We attempted to use a narrower slit, offset 
from the star by 0.6" to sample the wings of the star's p.s.f. This technique produced an 
unsaturated Ha profile, but it is probably too contaminated by surrounding ejecta to be 
trusted; see Subsection I 



Unfortunately r\ Car was not observed with Gemini/GMOS in 2005-2006 before 
the photometric fading. The ESO Science Archive does have spectra obtained with the 
VLT/UVES instrument in 200qj. We reduced the 2006 spectra with the standard UVES 



4 The 2006 UVES observations were obtained at ESO's Paranal Observatory under pro- 
gram ID 077.D-0618(A) by Weis et al. VLT/UVES observations of r) Car from 2002 to 2005 



are available a t http :/ /etacar.umn. edu/l The UVES spectral sequence is also described by 



Nielsen et al 
spectra. 



(120091 ) . but they did not note the potential problems with the later UVES 
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pipeline available from ESO. The spectra were extracted using a mesa function 3 by 2 pixels 
wide, about 0.75" by 0.5". The seeing was 0.8 to 0.9" so this extraction corresponds to 
about 1" on the sky. None of the 2006 spectra are concurrent with the decrease. The closest 
in time is from June 2006 but internal clues suggest that the slit was not centered on the 
star. The spectra do not closely resemble earlier UVES observations (2002 - 2004) which 
were definitely centered on the star; the stellar-wind features are weaker and the narrow 
emission lines are much weaker than what one expects to see in a ground-based spectrum. 
The published position and acquisition image for the June 2006 spectrum also suggest 
that the slit might be offset from the star. The June 2006 UVES spectrum thus appears 
unsuitable for our purposes here. The narrow emission lines are also weaker than expected 
in the two spectra from 2005 and the April 2006 spectrum and weaker than in the 2007 
GMOS spectra. This may also be due to pointing differences with respect to the earlier 
UVES spectra and the GMOS spectra. Although the April 2006 spectrum is referenced in 
the later discussion, it is used here with reservations for this reason. 



2.3. An unsolved puzzle: The effect of spatial resolution on spectroscopy of 

7] Car 

Provided the extraction size is less than 0.3", an HST spectrum of rj Car shows the 
stellar wind with only slight contamination. Ground-based spectroscopy, on the other hand, 
with seeing of the order of 1", also includes the spectrum reflected by dust at r ~ 0.15" 
to 1", plus narrow emission lines created in that region. Since reflection usually has little 
effect on the strengths of spectral features relative to the continuum, one might expect a 
ground-based spectrum of r] Car (i.e., central object plus reflected spectrum plus narrow-line 
emission) to closely resemble an HST spectrum, supplemented by t he narrow lines. In fac t, 



however, the earliest HST spectra showed that this is not the case ( [Davidson et al. 



19951 ) 
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Ground-based data show paradoxically weaker Balmer emission lines. 



Fig. [3] shows how conspicuous this effect is. Here we use VLT/UVES to exemplify a 
modern ground-based instrument, because r\ Car was observed with STIS and UVES on 
2003 February 12 and 14 respectively, only 1.2 days apart. We show the H<5 feature because 
it plays a leading role in §3.2. H5 emission in the UVES data has equivalent width 26 A 
compared to 33 A in the STIS spectrum (measured the same way in both data sets, with no 
correction for narrow lines). Old HST/FOS results and various ground-based data suggest 
that this difference is real, not merely an instrumental peculiarity of either STIS or UVES 



(IDavidson et al 



1995|). 



Why does the difference occur? Wavelength-dependent Balmer absorption within the 
reflecting dust-gas mixture seems implausible, because each stellar-wind emission line is 
much broader than the velocity dispersion in the relevant material. A more appealing 
explanation is that the spectrum "seen by" the dust intrinsically differs from what we see 
directly - in other words, the st ellar wind spectrum depends on the direction from which 
it is viewed 



Smith et al 



(120031 ) explored this idea via observations farther out in the 
Homunculus Nebula, and concluded that r\ Car's wind is strongly latitude-dependent. Note 
that the dusty 'Weigelt knots," which contribute strongly to any ground-based spectrum of 
the central object, are thought to be near the equatorial plane while our direct viewpoint is 
around 45° latitude e.g. Davidson et al (2001). 

This explanation is not entirely satisfying. It implies that the eq uivalent widths of 
Bal mer lines are substanti ally smaller at equatorial latitudes, whereas 
and 



Smith et al 



fl2003h 



Hillier fc Allenl (119921 ) suggested the same for polar directions. Apparently HST's 
high-spatial-resolution direct view of the stellar wind gives larger equivalent widths than 
we find in the reflected spectrum anywhere in the ejecta. (This statement probably applies 
to broad emission features in general, not just Balmer lines.) Furthermore, the UVES vs. 
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STIS comparison mentioned above seems to require a surprisingly large latitude variation 
in the H<5 equivalent width, almost a factor of two. Thus we fear that some other, more 
subtle process may be altering the reflected spectra. This problem merits further study. 

Independent of what causes the discrepancies, for this paper the practical implication 
is that one must be very careful when comparing ground- based spectra of rj Car to the 
HST/STIS data. In particular we cannot directly compare the equivalent widths. Moreover, 
the offset Ha observations mentioned in §2.2 are not useful here because they are most 
likely affected by reflection and emission in their spatial locales. In §4 below we compare 
recent Gemini/GMOS observations to earlier HST/STIS results, because neither data set 
spans two spectroscopic cycles; nor is any other existing data set suitable for this problem. 
In those comparisons, we exercise the necessary caution and caveats. 



3. Mid-Cycle Activity 

In recent years relatively little attention has been given to rj Car's behavior between 
spectroscopic events. At such times the companion star must be located 15-35 AU from 
the primary. Mid-cycle changes at UV-to-red wavelengths are therefore significant for at 
least one of two reasons: (1) Proposed explanations for the 5.54-year cycle generally do not 
predict irregular or short-term effects more than a few months from an event. If they do 
occur and are related to the cycle, then some new factor must be taken into account. (2) 
If, on the other hand, such effects are not related to the cycle, then they may represent 
phenomena in the primary wind. Observations in 1999-2002 and 2005-2007 therefore merit 
attention. Here we are concentrating on 2006-2007, supplemented by a few data from 
2000-2002 for comparison purposes. 

To avoid misunderstandings, note that continuous systematic changes occur in 
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mid-cycle, not directly related to the primary star. 



Mehneretal 



( 120 lOl ) have shown that 



some quasi-nebular spectral lines emitted close to our line of sight to the s tar gradually 



Damineli et al 



(120081 ) earlier found a 



brighten to a maximum near mid-cycle and then fade, 
similar result for one line of this type in an unresolved set of emission regions. The spectral 
features in question ([Ne III], [Ar III], [Fe III], etc.) probably indicate photoionization by 
the hot secondary star, and their behavior in the Mehner et al. analysis may represent 
the cyclic variation of column density between us and the orbiting secondary star. The 
most popular class of orbit-and-wind models can probably be adjusted to match this effect. 
Fluctuations described below, however, differ in several respects. They directly concern the 
spectrum of the primary stellar wind rather than nearby ejecta; they don't have a smooth 
obvious correlation with phase in the 5.5-year cycle; and each of them occurred in only one 
of two observed cycles. 



3.1. Photometry in 2006-2007 



The HST images confirm that the central star became fainter during 2006-2007 (Fig. 
[I]), but not as much as the ground-based photometry had led us to expect. From 1997 to 
2006 the star's apparent magnitude had brightened three times faster than the surrounding 
Homunculus, but its relative fading in 2006 was only about the same as for the Homunculus. 
Whatever the explanation is, this detail indicates that the change in 2006-2007 was not 
a simple reversal of the 1998-2006 trend. Moreover, the decline stopped and the upward 
brightening trend resumed after the 2009 spec troscopic event, leay i ng more or less intact 



our earlier speculations about secular changes (IMartin et al. 



2006b : 



Davidson et al. 



20051 ). 



The F330W filter showed a much steeper decline than F250W. This is interesting 
because the former samples Balmer continuum emission plus other continuum, whereas 
F250W includes many strong Fe II absorption features. During a spectroscopic event such 
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as 1998.0 or 2003.5, these absorptio n lines strengthen so much that the star becomes faint in 



the 230-290 nm wavelength range ( 


Cassatella. Giangrande. & Viotti 


1979: 


Altamore et al. 


1986; 


Viotti et al. 


1989; 


Davidson et al. 


1999b 


). This is why Fig. [TJ shows a substantial 



temporary decrease in F250W during the 2003.5 event, with scarcely any corresponding 
effect in F330W. Evidently the 2006-2007 fading was quite different. The same figure 
suggests two possibilities: (1) The increased F250W/F330W flux ratio in 2007 does not 
favor an explanation based on increased extinction by circumstellar dust. (2) Perhaps the 
Balmer continuum emission intrinsically weakened in 2006. 

The gr ound-based B, V, R, a nd I magnitudes all faded by about Am « +0.3 in 



2006-2007 flDiScala & Jones 



20081 ). While r\ Car's B and R values are strongly influenced 
by the very bright H/3 and Ha emission lines, V almost entirely measures continuum 
brightness. Unfortunately the central star could not be measured with HST/WFPC2 at 
these wavelengths, because every WFPC2 filter redward of 400 nm either saturates the 
central pixels in the shortest allowed exposure time, or else is affected by some particular 



emission line. 



3.2. Spectroscopy: 2007 vs. 2000-2002 

Motivated by the visual decline, its duration, and the peculiar UV photometry from 
HST/WFPC2, we obtained the Gemini/GMOS data described in Sj2Hto learn whether 
the spectrum had changed during the brightness decrease. Unfortunately no truly suitable 
earlier ground-based spectra are available for comparison - particularly none from 2005 or 
2006. Instead, therefore, we used STIS observations made in October 2001 and January 
2002, at approximately the same phase of rj Car's 5.54-year spectroscopic cycle as the 
GMOS data. There is no a priori reason to think that the 2006 photometric change was 
related to the cycle, but on the other hand it would be imprudent to use observations at 
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much different phases as comparisons. If, as usually assumed, the cycle is modulated by a 
hot companion star in an eccentric orbit, then Fig. [4] shows the approximate orbital location 
of the secondary star when In this part of the orbit the two 

stars are more than 20 AU apart and their motion is slow. 

The most conspicuous spectral features in 2007 (Fig. [2]) closely resembled the 2001-2002 
data. Consider for instance the brightest Balmer lines. Since Ha is overexposed in the 
GMOS data and H7 is confused with other features in this object, in Fig. [5] we show H/3 
and US. The measured equivalent widths for B.S turn out to be like those mentioned in 
§2.3: 26 A in the 2007 GMOS spectrum and 33 A in the 2001-2002 STIS data. Unlike the 
earlier examples, however, these values have been corrected for weak superimposed narrow 
lines. The internal uncertainties - chiefly systematic rather than statistical - are probably 
less than ± 2 A^ If we make allowances for the effect described in §2.3, evidently H£ was 



5 Most of the orbital parameters are extremely u ncertain, but Fig. H] is qualita ti vely valid 



for ne a rly all models that h ave been proposed. See 



Qkazaki et al. 



Humphreys fc Stanekl (120051 ); 



Ishibashi 



( 120081 ) and refs. therein. 



d2Q0JJ); 

6 Since the continuum of rj Car is hard to define and weaker additional emission lines 
contribute to the fluxes, these results depend on the measurement protocol. We measured 
"continuum" fluxes near 4080 and 4160 A and we integrated the net H£ flux from 4082 to 
4116 A (vac). The correction for unrelated emission on the long-wavelength wing is small, 
about —0.6 A in the GMOS data and less for STIS. In the last STIS data in early 2004 the 
equivalent width for B.S was only about 29 A, but lingering effects of the 2003.5 spectroscopic 



event may have reduced it below "normal" ( iDavidson et aL 



2005 ). Regarding e arlier epochs 



the first HST spectra gave an E.W. of 35 A for H<5 in 1991 ( IDavidson et al. 



19951 ) , or probably 



about 33 A if we could apply the corrections mentioned above. 
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essentially as strong in 2007 (relative to the continuum) as it had been 5.5 years earlier | 
The H/3 equivalent width was about 163, 175, and 156 A, respectively, in the 2001 and 2002 
STIS data and the 2007 GMOS data; this line may be less accurate than B.5 because its 
very large line/continuum flux ratio is disadvantageous for the measurements. 

But the H<5 profile does show one potentially significant difference: Weak P-Cygni 
absorption was present near —450 km s _1 in 2001-2002 and in the April 2006 UVES 
spectrum, but not in 2007. This is not merely a result of different spectral resolutions, 
because the STIS data in Fig. [5] have been smoothed to match the GMOS resolution. The 
difference in spatial resolution is probably not responsible either, since the STIS vs. UVES 
comparison in Fig. [3] shows at most a very weak P-Cygni f eature in the groun d-based data 



rather than STIS. (For reasons why that might occur, see ISmith et al.l (120031 ) along with 
§2.3 above.) The GMOS spectrum and B5 profile from February 2008 is essentially identical 
to the 2007 spectrum, while its July 2008 spectrum shows the reappearance of weak P Cygni 
absorption in B5 prior to the o nset of the 2009 event (Figure [6]). Although the spectroscopic 



analysis by iHillier et al. 



( 1200 ll ) pertained to the HST/STIS spectrum only one year after 
the 1998.0 event, they showed that the hydrogen emission originates predominantly in the 
outer wind of the primary star. Altogether, then, the P-Cygni difference is a real effect, 
meaning that the column density of excited hydrogen atoms along our line of sight was less 
after the photometric decline. 

A parallel, but more dramatic, difference occurred in some of the helium features. 
Helium lines whose lower levels are ls2s 1 S, ls2s 3 S, or ls2p 3 P° show prominent 
P-Cygni absorption in rj Car's mid-cycle spectrum; good examples are the triplet features 
at AA4026, 4472,4713. Qualitatively, one expects such absorption to occur in a dense 



7 Following the same protocol, the equivalent widths for B5 in the February and July 2008 
GMOS spectra are 25 A each and 26 A in the April 2006 UVES spectrum. 
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photoionized He + zone wher e recombination strongly populates the metastable ls2s levels 



( IQsterbrock fc Ferland 



20061 ). Such a zone is thought to exist in those parts of rj Car's 
primary wind that are closest to the hot second ary star before encountering the wind-wind 

(120081 ) and refs. therein. If Fig. H] is more or less 



shocks; see Section 6 of 



Humphreys et al. 



valid, that zone was situated between us and the primary wind in 2001-2002 and 2007; so 
P-Cygni absorption is unsurprising. For our purposes here, the interesting point is that the 
helium P-Cygni velocities differed between those two cycles. 

Fig. EH shows He I A4027. In 2001-2002 its P-Cygni feature was centered at V ~ -450 
km s _1 like that of H5, but in 2007 the absorption had broadened and shifted by about 
100 km s _1 to smaller (i.e., less negative) average velocities. This most likely indicates a 
quantitative change in the primary wind structure. Similar effects occurred in the other 
helium lines. The P Cygni absorption is also stronger in this line than the corresponding 
emission. We also note that the He I profile has continued to shift to even lower velocities 
in the February 2008 GMOS spectrum, but it has shifted bluewards, back to ~ —450 km 
s , in the July 2008 spectrum (Fig I7bl) . Interestingly, the He I A4027 P Cygni profile in 
the April 2006 UVES high resolution spectrum appears to have two absorption components 
centered at velocities of ~ —430 and ~ —400 km s -1 (Fig 17c!) . It is uncertain however, if 
the UVES slit is centered on the star or intercepting the same scattered and reflected light 
from the surrounding ejecta as the GMOS spectra. We plan to model the behavior of the 
He I lines during rj Car's spectroscopic cycle, but this is beyond the scope of this paper. 

By now it is clear that r\ Car is not a simple 5.54-year clock-like mechanism. HST/STIS 



data show a partial 



(IMartin et al. 



y different set of changes in the previous mid-cycle interval, 2000-2001 



2006d ). For example, P C ygni absorption in Ha, usually absent in a mid-cycle 



spectrum of the star (jSmith et al 



20031 ). showed a marked increase in absorption between 



March 2000, when there was no detectable P Cygni absorption, and October 2000. Its 
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appearance was brief since it was gone by January 2001. Unfortunately the October 2000 
STIS data only covered Ha so we could not compare the profiles of any of the other Balmer 
features or He I triplet lines. At the same time, there was a significant increase in the 
strength of the narrow -140 km s _1 absorption feature in Ha. The lower spectroscopic 
resolution of GMOS makes it difficult to determine if the narrow —140 km s _1 feature was 
present in 2007. 

The variable emission fea ture near 6307 A rea ches its maximum strength mid-cycle, in 



2001-2002 and again in 2005 flMartin et al 



2006c]). In the VLT/UVES spectra from 2005 



it had an equivalent width more than twice that previously observed. It is also present in 
the April 2006 UVES spectrum and June 2007 GMOS spectrum at 0.2 - 0.3A equivalent 
width_|, equal to its typical strength reported 1999-2002 between spectroscopic events. It is 
slightly stronger in the April 2006 spectrum relative to the June 2007 spectrum. We note 
that an increase in the UV flux from the central star occurs in 2001-2002 and 2005-2006 
and corresponds with the maximum strength of the 6307A line. 

There was at least one episode of mid-cycle activity in each of the last two spectroscopic 
cycles. Each episode was different and each occurred at a different phase relative to the 
spectroscopic events. In the next section we summarize the mid-cycle variability and discuss 
probable variations in the primary star's wind that would account for these changes. 



As outlined in ( Martin et al. 



2006d ) this line blends with an adjacent emission feature 
so it is difficult to accurately measure its strength without a more in depth analysis including 
spectra taken during a spectroscopic event. 
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4. Discussion 

The mid-cycle variability in 2006-2007 described in §3 occurred when the secondary 
star was at its greatest distance from the primary and most likely reveals basic fluctuations 
of rj Car's primary stellar wind, independent of the companion star and the 5.5-year cycle. 
This is the simplest view for reasons outlined below, and it makes sense in relation to several 
known or probable facts. For the first time in 160 years, the central star now strongly and 
directly affects ground-based photometry. The Homunculus nebula domina t ed the total 



Martin et al. 



fl2006bh : but 



apparent brightness throughout the 20th century, see refs. in 
since 1998 the star has been brightening rapidly and now accounts for half or more of 
the integrated light. This development may explain why ground-based magnitudes (i.e., 
Homunculus plus star) have recently been fluctuating with greater amplitude than in the 



1953-1995 record. See for example, the recent light curve at |http: / / etacar.fcaglp.unlp.edu.ar 



. Apparently our direct view of the star shows stronger variability than does the average of 
many viewpoints reflected in the Homunculus, though other factors should not be ignored. 

Since rj Car's spectroscopic events have attracted attention in recent years, there is 
a natural tendency to ascribe observed variability to the 5.5-year cycle, the wind-wind 
interaction, and the influence of the companion star. But in this case there are good 
arguments against that response. According to nearly all proposed models for the orbit, 
in 2006-2007 the wind-wind interaction should have been minimal because the companion 
star was then located far from the primary. During a 3-year interval around apastron, the 
separation is more than 20 AU (Fig. H]). This implies the following points for conventional 
models: 

1. Tidal effects are weaker than those at periastron by a factor of order 500. 

2. If the primary mass loss rate and line-of-sight wind speed are roughly 10~ 3 M Q yr _1 
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and wind speeds along the line of sight typically 500-600 km s l , then the optical 
depth for Thomson scattering at r > 20 AU is small, of the order of 0.1. 

3. Optical depths for continuum absorption are expected to be even smaller. A detailed 
model is beyond the scope of this paper, but this assertion is based on two generalities: 
Standard processes give r a b s < t sc at the likely densities and temperatures, and one can 
also judge the matter from the sizes of photospheres that produce suitable continua. 
The visual and near-UV photosphere in the primary wind cannot be very much larger 
than the circle for the primary in Fig. HJ R = 1 AU; otherwise its emergent continuum 
would be too cool. Therefore the optical depth for continuum absorption must be 



quite small at r ~ 20 AU. In conve ntional mode 



s the broad-emission- 



originates mostly within r < 6 AU (IHillier et al. 



2001 



Davidson et al- 



ine sp ectrum 



19951 ). where 



the secondary star and wind-wind interaction can have at most only a minor indirect 
influence. 



4. The radiation density is too high for appreciable c ircumstellar dust in this 



Likely grain formation distances are 150-300 AU (jDavidson &: Humphreys 



region. 



Kashi fc Soker 



19971). 



(]2008bl ) have suggested that dust can form close to the star, but also 
note that more than 170 days from periastron dust can form only at large distances. 
In any case as explained earlier dust, formation cannot explain the photometric 
decline. 

5. Along relevant photon paths (see point 6 below), optical depths through the fast, hot, 
low-density secondary wind are even smaller than those in the primary wind. 

6. The radiation density is too high for appreciable circumstellar dust in this regi on. 



Likely grain formation distances are 150-300 AU (jDavidson fc Humphreys! 1 199 71 ). 



7. When viewing Fig. HJ it is important to bear in mind the 45° inclination of the orbit 
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plane. At 2008.54, for instance, our line of sight to the primary star did not pass close 
to the secondary star even though the figure might give that impression. Between 
2005.5 and 2008.0 our projected view of the system would have changed slowly 
and continuously - at least in conventional models. The same statement applies to 
expected line-of-sight column densities. 

At visual and near-UV wavelengths, the shocked region between the two winds should 
be practically transparent for absorption and line formation. Its high temperatures, 
T > 10 6 K, greatly weaken free-free and other absorpti on opacities at thes e 



wavelengths. 



Okazaki et al 



igures depi cting shock simulati ons such as 



((20081), and 



Parkin et al. 



Pittard & Corcoran! (120021 ). 



(120091 ) do not depict the visual appearance; 



indeed the shocked regions are expected to emit less than 0.01% of r\ Car's total 
red-to- near-UV light. 

9. Flow simul ations show dramatic spiral patterns in the shocked str ucture around 



periastron (jOkazaki et al. 



2008 



Parkin et al. 



2009; 



Gull et al. 



20091 ). During most of 



the orbit, however, the secondary star moves at speeds less than 50 km s _1 , much 
slower than either wind; so the flow structure then approximates an axial symmetry 
around the the star-star axis, with only a modest tilt angle and negligible spirality. 

10. Near periastron the colliding wind struc ture tends to become unstable and chaotic 



(jMartin et al 



2006a 



Parkin et al. 



20091 ). but this effect is far less serious at larger 



distances and lower densities. 



11. 



he hot secondary star is thought to ionize the pre-shock zo nes in the primary wind 



(jDavidson 



1997 



Humphreys et al. 



2008 



Mehner et al. 



20101 ). This probably affects 



the He I P Cyg features described in §3, but it doesn't explain why the column 
densities and speeds should differ between two cycles. At most visual and near-UV 
wavelengths, continuum emission from this photoionized zone should be about two 
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orders of magnitude weaker than that of the primary wind. In the Balmer continuum 
it might account for almost 10% of the total, but this contribution should be quite 
steady, and insensitive to local gas densities etc. 

In summary, the above points plus the persistence of the spectroscopic changes supports 
variations in the primary's wind as opposed to a wind-wind interaction interpretation of the 
mid-cycle irregularities. We cannot absolutely rule out an explanation of that type, but the 
obvious alternative is far simpler - i.e., that the observations indicate fluctuations in the 
primary wind. Similar arguments apply to unconventional models wherein the secondary 
star plays the leading role, e.g. see Soker (2007) and Kashi and Soker (2009a). 

The decline in apparent brightness in 2006-2007 was temporary. The flux at all 
wavelengths subsequently leveled off and HST/WFPC2 photometry after the 2009 event 
(to be discussed in a forthcoming paper on the 2009 event) shows that the central star 
resumed its decades long increase in brightness. Furthermore the multi-wavelength visual 
photometry, and the near-UV flux ratio (Fig. [TJ show that increased extinction by dust 
cannot explain the decline. Allowable explanations for the observed decrease in brightness 
are thus limited. One feasible possibility is a decrease in the Balmer continuum emission 
which could be due a change in the stellar wind, assuming a simple dependence of the 
continuum on the wind density. 

The decrease in the H5 P Cygni absorption between 2001 and 2007 suggests a decrease 
in density of the primary star's wind between the two cycles at about the same phase in 
the orbit. If the density p(r) of an opaque wind decreases, then its photosphere tends 
to move inward. At constant luminosity (generally assumed for rj Car), this raises the 
characteristic temperature T c of the emergent radiation. Since T c is well above 10000° K, 
in this case, the brightness therefore tends to decrease at wavelengths greater than 300nm 
as the radiation shifts further in the UV. This is consistent with the observed photometric 
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decrease in 2006-2007. Crude estimates suggest that a density decrease on the order of 20% 
would suffice. The change in the He I profiles also show a corresponding decrease in the 
wind speed. The primary's wind thus became both slower and less dense mid-cycle. The 
Balmer and He I profiles show little change between the GMOS mid 2007 spectrum and the 
spectrum from early 2008. The primary's wind thus stayed in this slower, less dense state 
up until the time of the spectrum from mid 2008, just prior to the onset of the 2009 event. 



A possible explanation for these mid-cycle changes in the stellar w ind may be re 



r) Car's aspherical wind, confirmed with near-infrared interferometry (IWeigelt et al 



ated to 



20071 ). 



Our line of sight to rj Car is at an angle midway betwe en its fast, dense polar outflow and 



20031 ). A small shift in 



its slower, lower density wind in the equatorial region ( iSmith et al. 
7] Car's latitude-dependent wind could thus account for the observed spectroscopic changes 
with our line of sight intercepting more of the lower latitude wind from mid 2006 to mid 
2008. The He I double P Cygni absorption profile in the 2006 UVES spectrum is intriguing 
in this connection, but with doubts about the slit position, its interpretation is uncertain. 



Kashi fc Soker I (j2009bl ) have suggested that small changes in th e primary's wind can 



explain the short term variations observed in the X-ray light curve (ICorcoran 



20051 ) at 



apastron and propose that larger variations ma y have been respo nsible for the much shorter 



20091 ). However we note that 



X-ray minimum during the recent 2009 event (iPian et al. 
there were no large fluctuations in the X-ray flux during the time of the photometric decline 
and our observed spectral changes from w late 2006 to at least early 2008, although the 
obser ved 20% change in the wind velocity may result in only a 7% decrease in the X-ray 



flux ( jKashi fc Soker 



2009b|). 



We also noted mid-cycle changes in 2000-2002 although they were somewhat different 
and did not occur at the same phase in the orbit. These small mid-cycle fluctuations may 
be part of a longer term secular trend perhaps related to the star's rapid brightening in the 
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past decade as it recovers from the giant eruption. Most recent work on r] Car has focused 
on its spectroscopic events, but mid-cycle observations merit more attention. Mid-cycle 
variability can provide important information on the primary star's intrinsic instability that 
the spectroscopic events do not. 

We thank Jean-Rene Roy Deputy Director and Head of Science Operations of the 
Gemini Observatory, for granting us director's discretionary time on short notice to obtain 
spectra of rj Car in 2007. We also thank the staff and observers of the Gemini South 
observatory in La Serena for their help and support in putting together a successful 
observing program. We are also acknowledge use of the AAVSO International Database. 
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Table 1. Measured spatial sampling efficiency of 0.3" virtual aperture a for a point source 



HST camera F250W, F255W F330W, F336W 



ACS (HRC) b 0.594 ± 0.013 0.625 ± 0.010 
WFPC2 (PC) C 0.567 ± 0.030 0.669 ± 0.048 

a Weighting function l—r 2 /R 2 with R = 0.15". 



b Martin et al 



fl2006bh 



c The r.m.s. errors quoted here are relevant to com- 
parisons between ACS and WFPC2 data, but maybe 
pessimistic regarding brightness variations with a 
given filter. 
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Table 2. Photometry Results 









Exp Time Flux 






Dataset 


MJD 


Year 


(sec) Dcnsity c 


Magnitude a 


Average b 


ACS HRC/F250W Filter 


j9p603onq 


54120.426 


2007.055 


1.2 1.189 


6.212 


6.226±0.010 


j9p603opq 


54120.426 


2007.055 


1.2 1.159 


6.239 




j9p603orq 


54120.430 


2007.055 


1.2 1.174 


6.225 




j9p603otq 


54120.430 


2007.055 


1.2 1.172 


6.228 




WFPC2 HRC/F255W Filter 


ual40109m 


54335.188 


2007.643 


0.6 1.094 


6.302 


6.302 


ual40209m 


54510.797 


2008.124 


0.6 1.064 


6.333 


6.333 


ACS HRC/F330W Filter 


j9p603obq 


54120.410 


2007.055 


0.4 1.336 


6.085 


6.087±0.003 


j9p603odq 


54120.410 


2007.055 


0.4 1.338 


6.084 




j9p603ofq 


54120.414 


2007.055 


0.4 1.329 


6.091 




j9p603ohq 


54120.414 


2007.055 


0.4 1.334 


6.087 




WFPC2 HRC/F336W Filter 


ual40101m 


54335.121 


2007.643 


0.1 1.273 


6.138 


6.138±0.020 


ual40102m 


54335.125 


2007.643 


0.1 1.312 


6.105 




ual40103m 


54335.121 


2007.643 


0.1 1.253 


6.155 




ual40104m 


54335.125 


2007.643 


0.1 1.254 


6.154 




ual40201m 


54510.770 


2008.124 


0.1 1.154 


6.245 


6.262±0.019 


ual40202m 


54510.773 


2008.124 


0.1 1.135 


6.263 




ual40203m 


54510.773 


2008.124 


0.1 1.151 


6.247 




ual40204m 


54510.777 


2008.124 


0.1 1.105 


6.292 




a Magnitudc on the STMAG system. 








^The average 


and sigma 


of individual measurements in a set of 


exposures taken within a day 


of each other. 












c STMAG flux 


units are 


10 11 erg/cm 2 


/s/A. 
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Table 3. Spectral Observations 



Root Exp Length 



Name 


Year 


Grating 


(sec) 


Gemini South GMOS 


S20070616S0045 


2007.45 


B1200 


5.5 


S20070616S0046 


2007.45 


B1200 


5.5 


S20070618S0011 


2007.45 


B1200 


5.5 


S20070618S0012 


2007.45 


B1200 


5.5 


S20070630S0026 


2007.49 


R831 


3.5 


S20070630S0035 


2007.49 


R831 


1.5 


S20070630S0040 


2007.49 


R831 


4.5 


S20080211S0087 


2008.11 


B1200 


5.5 


S20080211S0088 


2008.11 


B1200 


5.5 


S20080718S0068 


2008.54 


B1200 


5.5 


S20080718S0069 


2008.54 


B1200 


5.5 


HST STIS 


o6ex03010 


2001.75 


G750M 


6.0 


o6ex03020 


2001.75 


G750M 


0.6 


o6ex03060 


2001.75 


G430M 


36.0 


o6ex030e0 


2001.75 


G430M 


36.0 


O 6ex030f0 


2001.75 


G430M 


3.0 


o(iex02010 


2002.05 


G750M 


6.0 


o6ex02020 


2002.05 


G750M 


0.4 


o6ex02030 


2002.05 


G430M 


18.0 


o6ex020u0 


2002.05 


G430M 


8.0 


o6ex020v() 


2002.05 


G430M 


36.0 



VLT UVES 



UVES. 2006-04-09 2006.27 CD#2 
T02:41:56.828 



60.0 
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Fig. 1. — Recent photometric behavior of rj Car. The curve marked "AAVSO" is a smoothed 
representation of groun d-based V-magnitudes for the Homunculus Nebula plus the star 



(IDiScala fc Jones II2008I ). Other data shown here include only the central star - or rather 
its opaque wind - observed with HST in the UV. Cro sses (+) show photo metry synthesized 



from STIS/CCD spectra near A ~ 330 and 25 nm (jMartin et al. 



measured in ACS/HRC images (jMartin et al 



2006bJ). Boxes (□) were 



(]2006bl ) and this paper). Asterisks (*) are 



photo metry from WFPC2 images. These HST points are STMAG values ( Koornneef et al. 



19861 ). An arrow marks the time when ground-based photometry indicated an unusual de- 
crease in brightness. The vertical dotted line indicates the time of the GMOS spectrum 
discussed in this paper. 



-34 - 




5000 5500 6000 6500 



1 
1 


u 


1 1 1 1 1 , 
i , , , , i 


i 
i 


i 

i 



5500 6000 6500 7000 7500 

WAVELENGTH IN VACUUM (Angstroms) 



Fig. 2. — A tracing of the entire GMOS spectrum of the star in July 2007. 
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Fig. 3.— The E5 profile in r] Car, observed with both HST/STIS and VLT/UVES at 
practically the same time in February 2003. The STIS tracing samples a diameter less than 
0.2", while the UVES spectrum represents a seeing-limited region about 1" across. Both 
data sets have been normalized so their continuum fluxes match near A « 4080 A and 4160 
A, and the bottom edge of the figure corresponds to zero flux. The spectral resolution of the 
UVES tracing has been degraded to match that of STIS, roughly 40 km s _1 . 
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V 



Fig. 4. — Orbital locations of the secondary star at times of the spectroscopy discussed here. 
The semi-major axis is about 17 AU and the eccentricity is assumed to be 0.85 (most authors 
favor the range 0.75-0.90), with periastron in early July 2003 and late January 2009. Our 
viewpoint is below the figure, except that the orbital plane is inclined by rough l y 45°. The 



apsidal orienta t ion s hown here is a compromise between estimates by 



Okazaki et al 



Ishibashil ( 120011 ) and 



( 120081 ). based on X-ray fluxes around apastron; it is uncertain by at least 



±20°. 
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Fig. 5. — Line profiles of two bright Balmer lines. The June 2007 GMOS spectrum is the 
solid black line. The solid gray and dashed black lines are the STIS spectra from October 
2001 and January 2002 respectively, smoothed to match the GMOS resolution of about 75 
km s _1 . 
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Fig. 6. — The H-delta line profiles for GMOS spectra observed July 2007 (solid black line), 
February 2008 (solid gray line), and July 2008 (dashed black line). Note that P-Cygni 
absorption is present in only the July 2008 spectrum. 
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Table 4. Instrument Parameters 





Slit Width 
(arcscconds) 


Slit Angle a 
(dog) 


Spectral 
R AA 


Gemini South 
GMOS 


0.50 


+ 160 


3744 b 
4396° 


HST STIS 


0.10 


+ 165 d 
-82 d 


153S0 b 
13400 c 


VLT UVES 


0.40 


+ 160 


80000 b 



a Thc slit angle is measured from north through east. All slits 
arc targeted on the central star. 

R for spectra blueward of 500 nm. 

C R for spectra rcdward of 500 nm. 

Slit orientation was +165 degrees in October 2001 and -82 
degrees in January 2002. 
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Fig. 7a. — A comparison of the velocity structure in the He I A4027 triplet line. As in Figure 
El the June 2007 GMOS spectrum is the solid black line. The solid gray and dashed black 
lines are the STIS spectra from October 2001 and January 2002 respectively smoothed to 
the same resolution as the GMOS spectrum. 
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Fig. 7b. — A comparison of the velocity structure in the He I A4027 triplet line in the 
GMOS spectra from July 2007 (solid black line), February 2008 (solid gray line), and July 
2008 (dashed black line). Note the blueward shift of the P-Cygni absorption in the July 
2008 spectrum. 
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Fig. 7c. — A comparison of the velocity structure in the He I A4027 triplet line in the GMOS 
spectra from July 2007 (solid black line) and the UVES spectrum from April 2006 (solid gray 
line). The UVES data has been binned by a factor of five to smooth out noise but it is still 
plotted with five times greater resolution than the GMOS spectra. Note that the P-Cygni 
absorption in the 2006 UVES spectra is asymmetric and has two velocity components. 



